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ABSTRACT 

The supernovae la data are used to analyze two general exact solutions for 
quintessence models. The best fit values for Q m o are smaller than in the A- 
term model, but still acceptable. With present-day data, it is not possible to 
discriminate among the various situations. 
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1. Introduction 

Recently, astronomers discovered an accelerated expansion of our Universe. It is well 
known that all known types of matter generate attraction, which leads to a decelerated 
expansion of the Universe. That discovery then reveals a new type of matter, which is now 
called quintessence or, sometimes, dark energy (Ostriker and Steinhardt 1995; Turner and 
White 1997; Chiba et al. 1997; Caldwell et al. 1998; Zlatev et al. 1999; Perlmutter et al. 
1999,b). 

The discovery of the presence of dark energy became possible when astronomers recog- 
nized that SN la can be the long expected standard candle for cosmological investigations. 
Two main features provide the use of SN la as a standard candle (Filippenko and Riess 

2000) : 

i) They are exceedingly luminous, comparable with luminosity of a whole galaxy; they 
can, thus, be detected and observed with high S/N ratio even at cosmological distances. 

ii) "Normal" SNe la have small variations among their peak absolute magnitudes (around 

0.3). 

The accelerated expansion of our Universe was discovered as a result of two projects: 
the High-z SN Search (Schmidt et al. 1998; Riess et al. 1998) and the Supernova Cosmology 
Project (Perlmutter et al. 1999, a). 

In fact, a new type of matter was predicted many years ago by A. Einstein, who included 
a A-term into his considerations (Einstein 1917). At the beginning of the past century, the 
A-term was just a new fundamental constant, and only much later it was really considered as 
a formidable challenge by both observational and theoretical cosmologists (Weinberg 1989; 
Dolgov et al. 1990; Zel'dovich 1992; Carroll et al. 1992; Carroll 2000; Sahni and Starobinsky 
2000; Rubakov 2000). Moreover, during the last 20 years cosmologists understood that 
this constant can be replaced with a scalar field, which induces the repulsive gravitational 
force dynamically. Accordingly, several models were proposed (Ratra and Peebles 1988; 
Peebles and Ratra 1988; Wetterich 1995, 1998; Copeland et al. 1998; Ferreira and Joyce 
1998; Liddle and Scherrer 1999; Steinhardt et al. 1999; Brax and Martin 1999; Sahni and 
Wang 2000; Binetruy 2000; Rubano and Scudellaro 2001), in order to explain the observed 
present acceleration of our Universe. 

Two of these models continued to be developed after the discovery of acceleration, and 
were also roughly elaborated and adapted for present-day data (Rubano and Scudellaro 

2001) . Here, we again use these models to such an end, but in a much more refined way: the 
goal is now to fit the observed data of apparent magnitude and redshift of the supernovae 
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Ia, and test the models themselves. 



2. Model description 

As said above, in this paper we discuss two models for quintessence, both based on a 
scalar field with a special type of potential. The field is minimally coupled with pressureless 
matter, and the total density parameter f2 of the Universe is fixed to be 1. A detailed 
discussion of the consequences of assuming such models in cosmology is given in (Rubano 
and Scudellaro 2001), so that we limit ourselves here only to a short summary of the results 
we need for our purpose. 

The main attractive feature of these models is that they allow a general exact solution 
of the field equations, obtained through a suitable transformation of variables. Anyway, 
independently of the fact that this is an exact solution, we also find that this solution 
reflects many properties of the real Universe correctly. 

The first model considers a potential of the form 

V{<p) = B 2 e^\ (1) 

where B 2 is a generic positive constant and a 2 is some fixed combination of universal con- 
stants 

»» = (2) 
c 

Actually, this kind of potential has already been widely discussed in the literature, but 
without any particular assumption on the value of a (Ratra and Peebles 1988; Peebles and 
Ratra 1988; Wetterich 1995, 1998; Copeland et al. 1998; Ferreira and Joyce 1998; Brax and 
Martin 1999; Sahni and Wang 2000; Binetruy 2000; Fabris et al. 2000). We stress that it 
is the particular choice of this constant given above that allows the exact integration of the 
field equations (see also (Barrow 1987; Burd and Barrow 1988)). 

The general solution of the cosmological equations (both for the metric and the scalar 
field) has five free parameters (including B 2 ) (Rubano and Scudellaro 2001). We fix two of 
them and keep three as free. But one of these three parameters just determines the present 
value of the scale factor of our Universe, which, in a spatially flat geometry, is not observable. 
It is not included into statistical analysis and does not affect the degrees of freedom of our 
analysis. 

We list below only the three cosmological functions which we need in our analysis (the 
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other ones can be found in (Rubano and Scudellaro 2001), of course) 

(l + z) =^^-y (3) 
~_ 2(1 + 2r 2 ) 

* " 3t.r(l + ^)' (4) 
1 + r 2 

fim = (1 + 2T 2 ) 2 ' (5) 

They are the redshift of the epoch, the Hubble parameter, and the fl m parameter of pres- 
sureless matter, and are expressed in terms of the dimensionless time r = t/t s . The free 
parameters are then the time scale t s and the present value of the dimensionless time tq. Lc 
us remark that t s is of the same order of magnitude (but not necessary equal to) as the age 
of the Universe. 

As to the second model, it considers a potential of the form 

V((p) = A 2 e^ + £V^, (6) 

with a 2 = 12Gtc/c 2 as before, and A 2 and B 2 free parameters. 

We have, now, one additional free parameter; therefore, according to the same consid- 
erations as above, we have to deal with three of them. 

The equations which describe the Hubble parameter, the redshift, and the density pa- 
rameter in the second model are 

s 3 A 2 sinh 2 r sin 2 r 

K l + z ) = , 2 ■ , 2 • 2 ' y') 

A sinn r-sm r 

^(sin(2r)-A 2 sinh(2r)) 
3(sm r — A smh r) 

_ 2(A 2 - l)(cos(2r) + A 2 cosh(2r) - 1 - A 2 ) 
mlj " (sin(2r)-A 2 sinh(2r)) 2 ' lJ 

The dimensionless time, in this case, is r = uot. Following Rubano and Scudellaro (2001), we 
use here u instead of t s , because of the fact that it is directly connected with the parameters 
in the potential of the scalar field, and has the meaning of a mass factor in theoretical 
considerations. So, the free parameters are u, A, and t . 
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In the analysis of the supernovae data, we use the bolometric distance. As explained 
better below, it can be expressed in terms of the "Hubble free" luminosity distance and of 
a parameter m , connected with the absolute magnitude and the Hubble parameter. The 
parameters of the first model (to and t s ) can be recasted into r and m . The parameters 
of the second model (A, To, and a;) can be recasted into A, To, and mo Once the best 
fit is made, it is easy to compute the relevant physical quantities H and fi m0 - In all the 
considerations below, H turns out to have the same value as in (Perlmutter et al. 1999, a). 
So, we concentrate on Q m o- 



3. SNe la Data 

The published data of the supernovae consist of 60 SNe la (Perlmutter et al. 1999, a). 
The data analysis and the determination of cosmological parameters can be considered in 
two steps. The first one is the measurement of the Hubble parameter for close supernovae 
(Calan - Tololo survey) (Hamuy et al. 1996), to be compared with the absolute magnitude 
M of a supernova SN la. The second step is the comparison of the high redshift supernovae 
with the theoretical prediction of bolometric distance: 

m = 5\og(D b ) + m ; (10) 

here, D b is the "Hubble free" bolometric distance 

and mo is a parameter connected to the absolute magnitude and the Hubble parameter. 

In data presented in (Perlmutter et al. 1999, a) there are several values for corrected ap- 
parent magnitude. Authors consider m p ^ ak and stretch luminosity corrected effective -B-band 
magnitude m e J . For the analysis of cosmological parameters only m e J is used, together 
with its errors a e //. 

m B 

There are several methods for SN la data analysis. Two of them are used in (Riess et 
al. 1998). The first one is the Multicolor Light Curve Shape (MLCS) method and the second 
one is a template fitting method. In (Perlmutter et al. 1999, a) another method is used. 

The data of both groups have the statistical errors approximately as a m ~ 0.25. 

We follow the authors of paper (Perlmutter et al. 1999, a) to analyze the models described 
in (Rubano and Scudellaro 2001). First of all, as a check of the procedure, we apply the 
flat cosmological model with a A-term to fit the data. The standard chi 2 algorithm of data 
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analysis reveals a good agreement of our analysis with published statistical values (Perlmutter 
et al. 1999, a). We use the complete set of data of 60 SNe la. It results x 2 = 1-75 per degree 
of freedom, not significantly different from x 2 — 1-76 found in (Perlmutter et al. 1999, a). 
The same is for the fl\ parameter. Since 4 points in the data are outliers, we can proceed 
with analysis and exclude these data from our considerations. The total number of SN la 
data then drops to 56. The x 2 P er degree of freedom in this case becomes 1.16, which is in 
good agreement with previously published results (Perlmutter et al. 1999, a) and is within 
la level. 



In our analysis we use the standard x 2 method. The analysis is done minimizing the 
value of weitghed x 2 '- 



where Wi is the weight of the i-th SN la, rrii is its 5-band effective apparent magnitude, 
and its magnitude as predicted with the models introduced before and thoroughly 

discussed in (Rubano and Scudellaro 2001). 



In the first model, it is possible to eliminate r from Eqs. (3) and (4), and to obtain an 
analytical expression for H(z). Thus, it is possible to compute m from Eqs. (10) and (11), 
and x 2 as a function of r and itlq. 

Firstly, we use 60 SN la data and get the x 2 minimum at m = 24.01, r = 1.04, with 
X 2 = 1-77 per degree of freedom. As it is unsatisfactory, we reject data which are out of the 
3a level, as done in (Perlmutter et al. 1999, a). 

After data rejection, the x 2 minimum drops down to x 2 — 1-195 per degree of freedom. 
It is definitely within one sigma level of the expected value of x 2 ■ The minimum now has 
other values than m = 23.985 and r = 1.268. 



If we accept the value of this minimum, we obtain, from Eqs. (4) and (5), Hq = 
7Qkms~ l Mpc~ l , n m0 = 0.15. 



4. Data analysis and fitting 




(12) 



4.1. 



The first model 



The situation is illustrated in Figs. 1-3. 



-7- 



4.2. The second model 

The second model has only been tested and fitted with 56 data of SNe la. The number 
of parameters in this case is equal to three. The true minimum of the x 2 is & t m o — 23.98, 
t = 0.8, and A = 1.182. We find a value of x 2 — 1-1906, which is definitely within one sigma 
level of expected value. After such value we obtain from Eqs. (8) and (9) that f2 m0 = 0.17. 

The x 2 value is a function of three arbitrary values: itlq, tq, and A. Therefore, the x 2 as 
a function of all parameters is impossible to plot, but we can nonetheless plot several slits. 

The situation is illustrated in Figs. 4-8. 

5. Conclusions 

In a quintessential universe we have analyzed the same data as in (Perlmutter et al. 
1999, a), where it is present only a cosmological constant, and found good values for x 2 m 
both cases. The values of Q m o found are rather different from the one found in (Perlmutter 
et al. 1999, a), but it is impossible to say if this is due to differences in the models or to 
influence of the measurement errors on the final values. 

In fact, in both models we have degeneracy in the parameters, particularly large in A 
(II model). This makes impossible to give significant confidence limits for the values of Q m o, 
which we found. Only very rough estimates can be given. Our main results are summarized 
in the following tables. 



Model I 



x 2 


m 


To 




t range 


Q m0 range 


1.195 


23.985 


1.268 


0.15 


0.82 1.40 


0.12 + 0.30 



Model II 



x 2 


m 


tq 


A 




Q m0 range 


1.19 


23.98 


0.8 


1.182 


0.17 


0.14 0.22 



As final remarks, we want to observe that our results are in a very good agreement with 
the one found in (Bahcall et al. 2000) in a completely independent way, and that the high 
degeneracy we get for the model parameters seems to support the opinion of those who claim 
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that it is very difficult to discriminate among theories on the basis of observational data only 
(Maor et al. 2000; Barger and Marfatia 2000). 
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Fig. 1. — The surface of \ 2 as a function of two variables, the To parameter and mo parameter. 
This function has a very definite and sharp minimum. 
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Fig. 2. — Countour plot of the levels of x 2 is shown. The one sigma level is the small 
white region on the graph. The two sigma level is the shadowed region. It reveals a large 
degeneracy of the x 2 function with respect to the parameter r . 
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Fig. 3. — The solid line corresponds to the slit of 2D \ 2 taken at the point of true minimum 
mo = 23.985. The dot-and-dashes curves corresponds to the slit taken at mo = 24.035, and 
the dashes curve to the slit taken at mo = 23.935. It results that the minimum in tq is 
strongly dependent on the m minimum. 




Fig. 4. — The surface of x 2 as a function of two parameters: mo and A; here, we fix To = 0.8. 
One can see the profile of this surface and minimum. 
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Tau0 = 0.8 




24.30 



Fig. 5. — Contour plot of the Fig. 4 surface. The la area is the small white region at 
m = 23.985 and at r = 1.268. The shadowed region is the 2a area, revealing that there is 
parameter degeneracy. 




Fig. 6. — We plot x 2 as a function of other two parameters, r and A, fixing m = 23.985. 
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Fig. 7. — Countour plot of Fig. 6 surface. The white region in the right corner below is the 
la level of parameters r and A. The shadowed region is the 2a level of these parameters. 
One can see again a large degeneration for them. 
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Fig. 8. — Three different slits at fixed m = 23.985 and r = 0.8, showing degeneracy in A. 



